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Abstract 

The excitation of H (Is — 2s) by electron impact in the presence and in the absence of 
the laser field is studied in the framework of the first Born approximation. The angular 
variation of the laser-assisted differential cross section (DCS) for atomic hydrogen by 
electronic impact is presented at various kinetic energies for the incident electron. The 
use of Darwin wave function as a semirelativistic state to represent the atomic hydrogen 
gives interesting results when the condition z/c <C 1 is fulfilled. A comparison with the 
non relativistic theory and experimental data gives good agreement. It was observed 
that beyond (2700 eV) which represents the limit between the two approaches, the 
non relativistic theory does not yield close agreement with our theory and that, over 
certain ranges of energy, it can be in error by several orders of magnitude. The sum rule 
given by Bunkin and Fedorov and by Kroll and Watson [23J has been verified in both 
nonrelativistic and relativistic regimes. 
PACS number(s): 34.80.Dp, 12.20.Ds 

1 Introduction 

Relativistic laser-atom physics has become recently as a new research area and constitute a 
new field of systematic experimental and theoretical study. In particular, the laser-asssited 
electron-atom scattering becomes a rapidly growing subject. These processes are however of 
fundamental interest and paramount importance, for instance in the laser heating of plasmas 
and high-power gaz lasers. One of its most remarkable features is the possibility of exciting 
the target via the absorption of one or more photons. In recent years considerable attention 
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has been given to these topics. Most of the works are theoretical and there is very little ex- 
perimental work at present. Several experiments have been performed, in which the exchange 
of one or more photons between the electron-target and the laser field has been observed in 
laser-assisted elastic [1] and inelastic scattering [2-5]. In particular, the excitation processes 
have been largely investigated in the literature by several authors [6-9], mainly in the pertur- 
bative weak-field limit. The first theoretical studies on the inelastic scattering were inspired 
from the pioneering works [10-12], in which the interaction between the free electron and the 
field can be treated exactly by using the exact Volkov waves [13]. Another investigation have 
been done for the exact analytical relativistic excitation 15*1/2 — > 1 5*1/2 of atomic hydrogen, 
by electron impact in the presence of a laser field, see [14]. For a summary of relativistic 
laser-atom collisions, see [15]. In the absence of the laser field, many authors have studied this 
process using numerical tools. Thus, Kisielius et al. [16] employed, the R-matrix method with 
nonrelativistic and relativistic approximations for the hydrogen like He + , Fe 25+ and U 91+ 
ions, where the case of transitions Is — > 2s and Is — > 2p as well as those between fine 
structure n = 2 levels was considered. Andersen et al. [17] have applied the semirelativistic 
Breit Pauli R-matrix to calculate the electron-impact excitation of the 2 5i/ 2 — > 2 P°/2 3/2 res " 
onance transitions in heavy alkali atoms. Payne et al. [18] have studied the electron-impact 
excitation of the 5s — > 5p resonance transition in rubidium by using a semi-relativistic Breit 
Pauli R-matrix with pseudo-states (close-coupling) approach. S. Taj et al. [19] have presented 
a theoretical semirelativistic model and have found that the semirelativistic Coulomb Born 
approximation (SRCBA) in a closed and exact form for the description of the ionization of 
atomic hydrogen by electron impact in the first Born approximation is valid for all geometries. 
The aim of this contribution is to add some new physical insights and to give the analytical 
formula for the excitation differential cross section for hydrogen from the Is ground state to 
the 2s excited state in the absence and in the presence of the laser field. Before we present the 
results of our investigation, we first begin by sketching the main steps of our treatment. For 
pedagogical purposes, we begin by the most basic results of our work using atomic units (a.u) 
in which one has (h — m e — e — 1), where m e is the electron mass at rest, and which will be 
used throughout this work. We will also work with the metric tensor g^ u = diag(l, —1,-1,-1) 
and the Lorentz scalar product which is defined by (a. b) = a M 6 M . The layout of this paper is 
as follows : the presentation of the necessary formalism of this work in section [2 and 3] , the 
result and discussion in section 4 and at last a brief conclusion in section 5. 
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2 Theory of the inelastic collision in the absence of de 
laser field 

In this section, we calculate the exact analytical expression of the semirelativistic unpolarized 
DCS for the relativistic excitation of atomic hydrogen by electron impact. The transition 
matrix element for the direct channel (exchange effects are neglected) is given by 



S fi = -i dt(ip p (x 1 )(pf(x2) | V d | ^ w (a?i)0i(ar 2 )) 



—i 



+00 



dt I rfr 1 ^ / (t,r 1 ) 7 V Pl (t,ri)(0 / (x 2 ) | V d | fafa)) (1) 



where 



V d = -~- (2) 

is the direct interaction potential, r x are the coordinates of the incident and scattered electron, 
r 2 the atomic electron coordinates, ryi = | r"i — r"2 | and r\ =\ ri |. The function ip Pi {xx) = 
ip p (t,ri) = u(p, s) exp(—ip.x)/\/2EV is the electron wave function, described by a free Dirac 
spinor normalized to the volume V, and 0jj(x 2 ) = <j>i,f{t, 1*2) are the semirelativistic wave 
functions of the hydrogen atom where the index i and / stand for the initial and final states 
respectively. The semirelativistic wave function of the atomic hydrogen is the Darwin wave 
function for bound states [17] . which is given by : 

0i(t,r 2 ) = exp(-i£ 6 (ls 1/2 )t)^i7(r 2 ) (3) 

where ^(lsi/ 2 ) is the binding energy of the ground state of atomic hydrogen and (pu( r 2) is 
given by : 

^S?(r 2 ) = (14 - ^«.V (2) )«( ± Vo(r 2 ) (4) 

it represents a quasi relativistic bound state wave function, accurate to first order in Z/c 
in the relativistic corrections (and normalized to the same order), with ip Q being the non- 
relativistic bound state hydrogenic function. The spinors u^' are such that uW = (1, 0, 0, 0) T 
and = (0, 1,0, 0) T and represent the basic four-component spinors for a particle at rest 
with spin-up and spin-down, respectively. The matrix differential operator a.V is given by : 



a.V 



( d Z d X ~idy\ 

d x + id y —d z 

d z d x — idy 

V d x + idy -d z 



(5) 
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For the spin up, we have : 



and for the spin down, we have 
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where 



N Dl = 2C/V4C 2 + 1 (8) 

is a normalization constant lower but very close to 1. Let us mention that the function r 2 ) 
in Eq. (1) is the Darwin wave function for bound states [18], which is also accurate to the order 
Z/c in the relativistic corrections. This is expressed as 0/(t,r 2 ) = exp(— iSb(2si/ 2 )t)Lp^\r 2 ) 
with £&(2si/ 2 ) as the binding energy of the 2sx/ 2 state of atomic hydrogen. 

/ 2 - r 2 ^ 
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for the spin up and 
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for the spin down. The transition matrix element in Eq. ([I]) becomes : 



Sfi 



dt / dr 1 dY 2 i) p (t, ri)7 ^(*, r i)0/(^, r 2 )(j)i(t, r 2 )V d 



and it is straightforward to get, for the transition amplitude, 



S 



fi 



2VjE f E i 



27rH inel (A)5(E f + S(2s 1/2 ) - E t - £{ls x/2 )) 



(12) 



where A = \p t — pf \ and 7 is the Dirac matrix. Using the standard technique of the QED, 
we find for the unpolarized DCS 



da 

dn~ f 



Ip/I 



\ Pi \ (4nc*y \2fr f 



\ u (Pf> s fh° u (Pi> s i)\ 2 \H inel {A)\' 



(13) 
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3 Calculation of the integral part 

The function Hi ne i(A) is found if one performs the various integrals : 

r+oo 

H inel (A) = / dv 1 e^I(r 1 ) (14) 



The quantity I{r\) is easily evaluated in the following way. We first write the explicit form of 

1 Z 



I n = / dr 2 4 s (r 2 ) 
Jo 



hs(r 2 ) (15) 



-ri2 n_ 

Next, we develop the quantity r^ 1 in spherical harmonics as 

1 . ^ YUfVY^ (r<)' 

Trr^l 2i + 1 < 16) 

where r> is the greater of r\ and r 2 , and r < the lesser of them. The angular coordinates of 
the vectors r x and r 2 are such that : fj = (0i,<^i) and r 2 = (02, ¥2)- We use the well known 
integral [T9] 

M m e" QU = -— re""* V ^4- ^e(a) > (17) 

then, after some analytic calculations, we get for I(ri) : 

,,6,1,414, 1 , 
^) = ^( ? -4) + ^ ? -- 5 (l + 8 ? )r 1 (18) 

The integration over ri gives rise to the following formula : 

Hi n ei(A) = | o + °° dr 1 e iAr U(r 1 ) = -^=(h + I 2 + J 3 ) (19) 

the angular integrals are performed by expanding the plane wave e* Ari in spherical harmonics 
as : 

e iA ri = ^4 7 rz^(Ar 1 )y, m (A)^(r 1 ) (20) 

/ m 

with A = pi — pj is the relativistic momentum transfer and A is the angular coordinates of 
the vector A. Then, after some analytic computations, we get for I\, J 2 and I3 the following 
result : 

4 r +oc , s„, . , . s 4 1 



4 r-r°° 3 
Jl = 27c 2 ,/ rfri r i e ~' ri io(Ar 1 ) 



27c 2 ((3/2) 2 + A 2 ) 

= ~ 4 » r ^ ^A(^) = 14 - *> ( ( 3/2)^ + A 2 ) 2 < 21 ' 

1 . A 2 -27/4 



((3/2) 2 + A 2 ) 3 
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It is clear that the situation is different than in elastic collision [2], since we have no singularity 
in the case (A — > 0) 

The calculation the spinorial part is reduced to the computation of traces of 7 matrices. 
This is routinely done using Reduce [20]. We consider the unpolarized DCS. Therefore, the 
various polarization states have the same probability and the actual calculated spinorial part 
is given by summing over the final polarization Sf and averaging aver the initial polariztion 
Si. Therfore, the spinorial part is given by : 

^Kp /)S/ ) 7 V^, S ,)| 2 = Tr{ 7 °(Ac + c 2 ) 7 %c + c 2 )} 

s i s f 

= 2c 2 [^-( W ) + c 2 ] (22) 

We must, of course, recover the result in the nonrelativistic limit (7 — > 1), situation of which 
the differential cross section can simply given by : 



da 


K/ 
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dVL f 
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with |A nr | = |Kj — Kj| is the nonrelativistic momentum transfer and the momentum vectors 
(Kj, K.j) are related by the following formula : 

K^dK^-S/^ 1 / 2 (24) 



4 Theory of the inelastic collision in the presence of de 
laser field 



The second-order Dirac equation for the elctron in the presence of an external electromagnetic 
field is given by : 



V -lV r J -/• 



2c 



ip(x) = 



(25) 



where = ^[7^,7^] is the tensors related to Dirac matrices 7^ and F^ v = d^Ay — d v A^ i 



is 



the electromagnetic field tensor. N 1 is the four-vector potential. The plane wave solution of 
the second-order equation is known as the Volkov state [13J. 



( 1 + J4J\ exp 



-i(qx) — i 



c(kp) 



(26) 
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We turn now to the calculation of the laser-assisted transition amplitude. The instantaneous 
interaction potential is given by 

V d = — -- (27) 
ri2 ri 

where ri 2 = \t± — T2\; ri are the electron coordinates and T2 are the atomic coordinates. 
The transition matrix element corresponding to the process of laser assisted electron-atomic 
hydrogen from the initial state % to the final state / is given by 



Sfi = -ij dt (Vv(ri)0/(r 2 )| WJrO^rs)) (28) 

Proceeding along the lines of standard calculations in QED, one has for the DCS 

da + °° r/7f (s) 

£ ?r ( 29 ) 



d ^f st^oo dttf 

with 

da {s) _ |q/| 1 / 1 ,„,(.) ,2^ itt /a \ i - 



<Kl f | qu| (4vrc 2 ) 2 I 2 



J2\M^\ 2 \\H mel (A s 



(30) 

Qf=Qi+SU+E ls y 2 -E 2sl /2 



The novelty in the various stages of the calculations is contained in the spinorial part ~ J2 Si J2 Sf \M' 
that contains all the information about the spin and the laser-interaction effects. This quan- 
tity can be obtained using REDUCE (20] and are explicitly given in our previous work [2T] . 
After some analytical calculations, the integral part H ine i(A s ) reduces to 

H inel (A s ) = -%[h{s) + h(s) + I 3 (s)] (31) 



n 12 



with 



his) 



1 



27c 2 ((3/2) 2 + A s ' 



2-i 



/2(s) = h^-^m^f- (32) 



9 V 8c 2; ((3/2) 2 + A s 2 ) 3 
with A s = qj — qj — sk is the momentum transfer in the presence of the laser field. 



8 







Probability density for radial Darwin (1s) 

Probability density for radial Dirac 1s,, 2 



5x1 J " 




Atomic number Z 

Figure 1: Behavior of the probability density for radial Darwin wave function compared with 
that of the Dirac wave function for small distances and for increasing values of the atomic 
charge number. 

5 Results and discussions 
5.1 In the absence of the laser field 

In presenting our results it is convenient to consider separately those corresponding to non- 
relativistic regime (the relativistic parameter 7 ~ 1) and those related to relativistic one (the 
relativistic parameter 7 ~ 2). Before beginning the discussion of the obtained results, it is 
worthwhile to recall the meaning of some abbreviation that will appear throughout this section. 
The NRDCS stands for the nonrelativistic differential cross section, where nonrelativistic plane 
wave are used to describe the incident and scattered electrons. The SRDCS stands for the 
semirelativistic differential cross section. 

We begin our numerical work with the study of the dependence of the probability density 
for radial Darwin and Dirac wave functions on the atomic charge number Z . 

So long as the condition Za <C 1 is verified, the use of Darwin wave function do not have 
any influence at all on the results at least in the first order of perturbation theory. So, the 
semirelativistic treatment when Z increases may generate large errors but not in the case of 
this work. In this paper, we can not have numerical instabilities since there are none. For the 
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Scattering angle e (degree) Scattering angle e (degree) 



Figure 2: The long-dashed line represents the 
semi-relativistic DCS, the solid line represents 
the corresponding non-relativistic DCS for a 
relativistic parameter (7 = 1.5) as functions of 
the scattering angle 9. 



Figure 3: The solid line represents the semi- 
relativistic DCS, the long-dashed line repre- 
sents the corresponding non-relativistic DCS 
for various values of the relativistic parameter 
(7 = 1.5, 7 = 2 and 7 = 2.5) as functions of 
the scattering angle 9. 



sake of illustration, we give in figure 1 the behavior of the probability density for radial Darwin 
wave functions as well as that of the exact relativistic Dirac wave functions for different values 
of Z. As you may see, even if it is not noticeable in figure 1, there are growing discrepancies 
for Z = 10 and these become more pronounced when Z = 20. The QED formulation shows 
that there are relativistic and spin effects at the relativistic domain and the non relativistic 
formulation is no longer valid. 

In the relativistic regime, the semirelativistic differential cross section results obtained 
for the Is — > 2s transition in atomic hydrogen by electron impact, are displayed in figures 2 
and 3. In this regime, there are no theoretical models and experimental data for comparison 
as in nonrelativistic regime. In such a situation, it appears from figures [2] and [3] that in the 
limit of high electron kinetic energy, the effects of the additional spin terms and the relativity 
begin to be noticeable and that the non-relativistic formalism is no longer applicable. Also a 
peak in the vicinity of Of = 0° is clearly observed. 

The investigation in the nonrelativistic regime was carried out with 7 as a relativistic pa- 
rameter and 9 as a scattering angle. In atomic units, the kinetic energy is related to 7 by 
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Scattering angle e (degree) Scattering angle e (degree) 



Figure 4: The solid line represents the semi- Figure 5: The variation of the SRDCSs with 

relativistic DCS, the long-dashed line repre- respect to 6, for various kinetic energies, 

sents the corresponding non-relativistic DCS 
for a relativistic parameter 7 = 1.00053 as 
functions of the scattering angle 9. 



the following relation : J5& = c 2 (7 — 1). Figure 4 shows the dependence of DCS, obtained in 
two models (SRDCS, NRDCS), on scattering angle 9. In this regime, it appears clearly that 
there is no difference between these models. Figure 5 shows the variation of the SRDCS with 
9 for various energies. It also shows approximatively in the interval [-5, 5] that the SRDCS 
increases with 7, but decreases elsewhere. Figure 6 presents the observed and calculated an- 
gular dependence of Is — 2s differential cross section of e~ — H scattering at incident energie 
200 eV. Results obtained in the approaches (semirelativistic and nonrelativistic approxima- 
tions) are indistinguishable and in good agreement with the experimental data provided by J. 
F. Williams [22]. 

5.2 In the presence of the laser field 

For the description of the scattering geometry, we work in a coordinate system in which k\\e z . 
This means that the direction of the laser propagation is along Oz axis. We begin by denning 
our scattering geometry. The undressed angular coordinates of the incoming electron are 
9i = 90°, (pi = 45° and for the scattered electron, we have 0° < 9 f < 180° and (fi f = 45°. 
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Figure 6: The variation of the differential Is — 2s cross section of e~ — H scattering at 200 eV. 
The dots are the observed values of J. F. Wiliams (1981) ; the solid line represents the semi- 
relativistic approximation and the long-dashed line corresponds to the non-relativistic DCS. 

When we take the zero electric field strength (E = a.u), one finds overlapping curves for 
the tree approaches [SRDCS (without laser), SRDCS (with laser) and the Non relativistic 
DCS]. It represents our first consistency check of our calculations and it is shown in figure [71 
Figure |8] shows the envelope of the laser assisted SRDCS as a function of the net number of 
photons exchanged. The cutoffs are s ~ —50 photons for the negative part of the envelope 
and s ~ +50 photons for the positive part. In figure El for the geometry 9i = 90°, fa = 45° 
and (ftf = 45°, we have made simulations concerning the laser-assisted SRDCS for a set of net 
number photons exchanged. These sets (±10, ±20, ±30, ±40, ±45 and ±50) show that at 
±50 the SRDCS (with laser) is almost close to the SRDCS (without laser). But for instance 
at ±10, we have several orders of magnitude as a result of the difference between the two 
approaches. We return to the first case ±50, the convergence reached here is called the sum- 
rule that was shown by Bunkin and Fedorov as well as by Kroll and Watson [23J. The figures 
flH] and |9]) correctly introduce correlation in the net number of photons exchanged that reaches 
the well-known sum-rule. As you can see from the figure [H the SRDCS full off abruptly beyond 
the interval [—50, +50] and figure [9] shows clearly that beyond ±50, the sum-rule is clearly 
checked. Figure [IU] show clearly that for the relativistic regime (7 = 2.0 and S = 0.05 a.u) 
and for the geometry (6>j = 45°, fa = 0°, <fif = 45° and 9f — 45°), the value of the cutoffs 
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Figure 7: The various DCSs as a function of Of for an electrical field strength of £ = a.u, a 
relativistic parameter 7 = 1.0053 




Figure 8: Envelope of the SRDCS in the 
non relativistic regime (7 = 1.0053 and £ = 
0.05 a.u) for a number of net photons ex- 
changed ±100. 



Figure 9: The variation of the SRDCSs with 
and without laser for various numbers of pho- 
tons exchanged (±10, ±20, ±20, ±40, ±45) 
and ±50). 
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Figure 10: Envelope of the SRDCS in the relativistic regime (7 = 2.0 and 8 = 0.05 a.u) for a 
number of net photons exchanged ±30000. 

have been changed to s — —13986 for the negative part of the envelope and s = +13912 
for the positive part. In this regime, when the number of photons exchanged increases, the 
convergence of the two SRDCSs (with and witout laser) will be reached at approximatively 
±13945, but as our computational capacity is limited, it is not possible, actually, to give the 
figure which illustrates such a situation. 

6 Conclusion 

This paper display the results of a semirelativistic excitation of atomic hydrogen by electronic 
impact in the presence and absence of the laser field. The simple semirelativistic Darwin 
wave function that allows to obtain analytical results in an exact and closed form within the 
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framework of the first Born approximation is used. Our results have been compared with 
previous nonrelativistic results, revealing that the agreement between the different theoretical 
approaches is good in the nonrelativistic regime. The nonrelativistic treatment has been shown 
to be no longer reliable for high energies. It has been also found that a simple formal analogy 
links the analytical expressions of the unpolarized differential cross section in the absence of 
the laser field and the laser-assisted unpolarized differential cross section. 
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